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1. Introduction 

Complex, self- regulating systems such as the human 
heart must process inputs with a broad range of charac- 
teristics to change physiological data and time series. 13 
Many of these physiological time series seem to be highly 
chaotic, represent nonstationary data, and fluctuate in 
an irregular and complex manner. One hypothesis is that 
the seemingly chaotic structure of physiological time se- 
ries arises from external and intrinsic perturbations that 
push the system away from a homeostatic set point. An 
alternative hypothesis is that the fluctuations are due, 
at least in part, to the systems underlying dynamics. In 
this review, we describe new computational approaches 
based on new theoretical conceptsfor analyzing physio- 
logical time series. Well show that the application of 
these methods could potentially lead to a novel diagnos- 
tic tool for distinguishing healthy individuals from those 
with congestive heart failure (CHF). 

Physiological Time Series 

Recent research suggests that physiological time series 
can possess fractal and self-similar properties, which are 
characterized by the existence of long-range correlations 
(with the correlation function being a power-law type). 
However, until recently, the analysis of such fluctuations 
fractal properties was restricted to computing certain 
characteristics based on the second moment of the data, 
such as the power spectrum and the two-point autocorre- 
lation function. These analyses indicated that the fractal 
behavior of healthy, free-running physiological systems 
could be characterized, at least in some cases, by 1/f-likc 
scaling of the power spectra over a wide range of time 
scales. 4,5 A time series that exhibits such long-range cor- 
relations with a power-law correlation function and is also 
homogeneous (different parts of the scries have identical 
statistical properties) is called a monofractal series. How- 
ever, many physiological time series are inhomogeneous 
in the sense that distinct statistical and scaling proper- 
ties characterize different parts of the series. In addition, 
there is some evidence that physiological dynamics can 
exhibit nonlinear properties. 6-12 Such features are of- 



ten associated with multifractal behaviorthc presence of 
long-range power-law correlations in the higher moments 
of the time serieswhich, unlike monofractals, are nonlin- 
ear functions of the second moments scaling exponents. 13 
Up until recently, though, robust demonstration of mul- 
tifractality of nonstationary time series was hampered 
by problems related to significant bias in the estimates 
of the datas singularity spectrum, due to the time scries 
diverging negative moments. A new wavelet-based multi- 
fractal formalism helps address such problems. 13 Among 
physiological time series, the study of the statistical prop- 
erties of heartbeat interval sequences has attracted much 
attention. 14-17 The interest is at least partly due to the 
facts that the heart rate is under direct neuroautonomic 
control; interbeat interval variability is readily measured 
by noninvasive means; and analysis of heart-rate dynam- 
ics could provide important diagnostic and prognostic 
information. Thus, extensive analysis of interbeat in- 
terval variability represents an important quantity for 
elucidating possibly nonhomeostatic physiological vari- 
ability. Figure 1 shows examples of cardiac interbeat 
time series (the output of a spatially and temporally in- 
tegrated neuroautonomic control system) for healthy in- 
dividuals and those with CHF. In the conventional ap- 
proaches to analyzing such data, we would assume the ap- 
parent noise has no meaningful structure, so we wouldnt 
expect to gain any understanding of the underlying sys- 
tem through the study of such fluctuations. Conventional 
studies that focus on averaged quantities therefore usu- 
ally ignore these fluctuations in fact, theyre often labeled 
as noise to distinguish them from the true time series 
of interest. However, by adapting and extending meth- 
ods developed in modern statistical physics and nonlin- 
ear dynamics, the physiological fluctuations in Figure 1 
can be shown to exhibit an unexpected hidden scaling 
structure. 5,10,13 ' 18,19 Moreover, the fluctuation dynamics 
and associated scaling features can change with patho- 
logical perturbation. These discoveries have raised the 
possibility that understanding the origin of such tempo- 
ral structures and their alterations through disease could 



elucidate certain basic aspects of heart-rate control mech- 
anisms and increase the potential for clinical monitoring. 
But despite this considerable progress, several interest- 
ing features must still be analyzed and interpreted. The 
theoretical concepts we discuss here are based on the pos- 
sible Markov properties of the time series; a cascade of 
information from large time scales to small ones that are 
built based on increments in the time series; and the 
extended self-similar properties of the beat-to-beat fluc- 
tuations of healthy subjects as well as those with CHF. 
The method we describe uses a set of data for a given 
phenomenon that contains a degree of stochasticity and 
numerically constructs a relatively simple equation that 
governs the phenomenon. In addition to being accu- 
rate, this method is quite general, can provide a ratio- 
nal explanation for complex features of the phenomenon 
under study, and requires no scaling feature or assump- 
tion. As we analyze cardiac interbeat intervals, well also 
look at new methods for computing the Kramers-Moyal 
(KM) coefficients for the increments of interbeat intervals 
fluctuations, Ax(t) = [x(t + r) — x(t)]/a T , where <r T is the 
standard deviations of the increments in the interbcats 
data. Whereas the first and second KM coefficients (rep- 
resenting the drift and diffusion coefficients in a Fokkcr- 
Planck [FP] equation) have well-defined values, the third- 
and fourth-order KM coeffi- cients might be small. If so, 
we can numerically construct an FP evolution equation 
for the probability density function, P(Ax,t) which, in 
turn, can be used to gain information on the PDFs evo- 
lution as a function of the time scale r. 2023 

Regeneration of Stochastic Processes 

Lets start by examining the computations that lead 
to the numerical construction of a stochastic equation. 
This equation describes the phenomenon that generates 
the data set, which is then utilized to reconstruct the 
original time series. Two basic steps are involved in the 
numerical analysis of the data and their reconstruction. 
Data Examination We must first examine the data to 
see whether they follow a Markov chain and, if so, we 
estimate the Markov time scale tw As is well known, 
a given process with a degree of stochasticity can have 
a finite or an infinite Markov time scale, which is the 
minimum time interval over which the data can be con- 
sidered to be a Markov process. 20,24 To determine the 
Markov scale tM, we note that a complete characteriza- 
tion of the statistical properties of stochastic fluctuations 
of a quantity x{t) requires the numerical evaluation of 
the joint PDF P n (x\, t\\ • ■ • ; x n , t n ) for an arbitrary n, 
the number of the data points in the time series x(t). If 
the time series x(t) is a Markov process, an important 
simplification can be made as P n , the n-point joint PDF, 
is generated by the product of the conditional probabili- 
ties P{xi+\, U + i \xi, ti), for i = 1, • ■ • , n — 1. A necessary 
condition for x(t) to be a Markov process is that the 




FIG. 1. Interbeats fluctuations of healthy subjects (top), 
and those with congestive heart failure (bottom). 

Chapman-Kolmogorov (CK) equation, 

P(x 2l t 2 \xi,t 1 ) = J d(X3) P{x2,h\X3,t3) P(x3,tz\xi,tl) , 

(1) 

should hold for any value of t 3 in the interval t 2 < 
t 3 < t\. One should check the validity of the CK equa- 
tion for various x\ by comparing the directly-computed 
conditional probability distributions P(x 2 , t 2 \xi, t\) with 
the ones computed according to right side of Eq. 
(1). The simplest way to determine tm for sta- 
tionary or homogeneous data is the numerical com- 
putation of the quantity, S = \P{x 2 ,t 2 \x\,ti) — 
J dx 3 P(x 2l t 2 \x 3l t 3 ) P(x 3l t 3 \xi,ti)\, for given xi and x 2 , 
in terms of, for example, t 3 — t\ (taking into account 
the possible numerical errors in estimating S). Then, 
tM =t$—ti for that value of t 3 — 1\ for which S vanishes 
or is nearly zero (achieves a minimum). (2) Numeri- 
cal construction of an effective stochastic equation that 
describes the fluctuations of the quantity x{t), represent- 
ing the time series, constitutes the second step. The CK 
equation yields an evolution equation for the PDF P(x, t) 
across the scales t. The CK equation, when formulated 
in differential form, yields a master equation which takes 
the form of a FP equation: 
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The drift and diffusion coefficients, D^(x,t) and 
D^ 2 \x,t), are computed directly from the data and the 
moments of the conditional probability distribu- 

tions: 
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FIG. 2. Test of Chapman-Kolmogorov equation for 
xi = —5, Xi = and x\ — 5. The bold and open sym- 
bols represent, respectively, the directly-evaluated PDF and 
the integrated PDF. The PDFs are shifted in the vertical di- 
rections for better presentation. Values of x are measured in 
units of the standard deviation. 

Note that the above FP formulation is equivalent to the 
following Langevin equation: 

d 



dt 



x(t)=D^(x) + ^DW(x) /(f) 



(5) 



where f(t) is a random force with zero mean and Gaus- 
sian statistics, ^-correlated in i, i.e., (f(t)f(t')) = 2S(t — 
t 1 ). We note that the numerical reconstruction of a 
stochastic process does not imply that the data do not 
contain any correlations, or that the above formulation 
ignores the correlations. 

Equation (5) enables us to reconstruct a stochastic time 
series x(t) 7 which is similar to the original one in the sta- 
tistical sense. The stochastic process x(t) is regenerated 
by iterating Eq. (5) which yields a scries of data with- 
out memory. To compare the regenerated scries with the 
original x(t), we must take the temporal interval in the 
numerical discretization of Eq. (5) to be unity (or renor- 
malize it to unity, if need be). However, the Markov time 
is typically greater than unity. Therefore, we correlate 
the data over the Markov time scale tM, for which there 
are a number of methods. 21 ' 22,25 A new technique that 
we have used in our own studies, which we refer to as the 
kernel method, is one according to which one considers a 
kernel function K{u) that satisfies the condition that, 



K{u)du = 1 , 



(0) 



such that the data are determined, or reconstructed, by 
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where h is the window width. For example, one of the 
most accurate kernels is the standard normal density 




FIG. 3. The drift and diffusion coefficients (x) and 
D^ 2 '(x), estimated by Eq. (3). For the healthy subjects 
(triangles) D^(x) and D^ 2 \x) follow linear and quadratic 
behavior in x, while for patients with CHF (squares) they 
follow third- and fourth-order behavior in x. 



function, K(u) = (27r) _1 / 2 exp(— \u 2 ). In essence, the 
kernel method represents the the time series as a sum 
of "bumps" placed at the "observation" points, with its 
kernel determining the shape of the bumps, and its win- 
dow width h fixing their width. It is evident that, over 
the scale h, the kernel method correlates the data. 

3. Analysis of Fluctuations in Human Heartbeats 

To show how the above reconstruction method is used 
in practice, and to demonstrate its utility, the method 
has been applied to reconstruction of the fluctuations in 
human heartbeats of both healthy and ill subjects by 
taking h ~ tM- Recent studies 5 ' 19 ' 26 ' 27 reveal that under 
normal conditions, beat-to-beat fluctuations in the heart 
rates might display extended correlations of the type typ- 
ically exhibited by dynamical systems far from equilib- 
rium. It has been argued, 26 for example, that the various 
stages of sleep might be characterized by extended cor- 
relations of heart rates separated by a large number of 




FIG. 4. The curves show, from top to bottom, the actual 
interbeat data (for a healthy subject), the regenerated data 
using the corresponding Langevin equation, and the regener- 
ated data using the kernel method. The time series are shifted 
in the vertical directions for better presentation. 



beats. While the existence of extended correlations in the 
fluctuations of human heartbeats is an interesting and 
important result, we argue that the Markov time scale 
tM and the associated drift and diffusion coefficients of 
the interbeat fluctuations of healthy subjects and those 
with CHF help one to better distinguish the two classes of 
subjects, particularly at the early stages of the disease, 
as these quantities have completely different behaviour 
for the two classes of subjects. 

Both daytime (12:00 pm to 18:00 pm) and nighttime 
(12:00 am to 6:00 am) heartbeat time series of healthy 
subjects, and the daytime records of patients with CHF 
have been analyzed by this method. 45,46 The data base 
includes 10 healthy subjects (7 females and 3 males with 
ages between 20 and 50, and average age of 34.3 years), 
and 12 subjects with CHF (3 females and 9 males with 
ages between 22 and 71, and average age of 60.8 years). 
Figure 1 presents the typical data. 

As the first step, the Markov time scale tu of the data 
is computed. From the daytime data for healthy subjects 
the values of tu are computed to be (all the values are 
measured in units of the average time scale for the beat- 
to-beat times of each subject), t.M = 3, 3, 3, 1, 2, 3, 3, 2, 3 
and 2. The corresponding values for the nighttime 
records are, tM arc 3,3,1,3,3,2,3,3,2 and 3, respec- 
tively, comparable to those for the daytime. On the 
other hand, for the daytime records of the patients 
with CHF, the computed Markov time scales are, t^j — 
151,258,760,542,231,257,864,8,366,393,385, and 276. 
Therefore, the healthy subjects are characterized by tM 
values that are much smaller than those of the patients 
with CHF. Thus, one has an unambiguous quantity for 
distinguishing the two classes of patients. 

Next, the validity of the CK equation for describing the 
phenomenon is checked for several Xi-triplets by compar- 



ing the directly-computed conditional probability distri- 
butions P{x2,t2\x\,t\) with the ones computed accord- 
ing to right side of Eq. (1). Here, x is the interbeat and 
for all the samples we define, x = (x — x)/a, where x and 
a are the mean and standard deviations of the interbeat 
data. In Figure 2, the two PDFs, computed by the two 
methods, are compared. Assuming the statistical errors 
to be the square root of the number of events in each bin, 
the two PDFs are statistically identical. 

The corresponding drift and diffusion coefficients, 
D^\x) and D^(x), are displayed in Figure 3, demon- 
strating that, in addition to the Markov time scale 
t]\j, the two coefficients provide another important in- 
dicator for distinguishing the ill from healthy subjects: 
For the healthy subjects the drift Z)W and the diffu- 
sion coefficients D^ 2 \x) follow (approximately) linear 
and quadratic functions of x, respectively, whereas the 
corresponding coefficients for patients with CHF follow 
(approximately) third- and fourth-order equations in x. 
Thus, for the healthy subjects, 



-0.12a; 



(8) 



D {2) (x) = (5 - 4.2a: + 7x 2 ) x 10~ 2 , (9) 
whereas for the patients with CHF, 

D {1) {x) = -(26a- + 18a- 2 + 7a; 3 ) x 10 -4 , (10) 



D (2) (x) = (6- 7x + 5a- 2 + 3a; 3 + 2a; 4 ) x 10' 



(11) 



which were also obtained by Kuusela. 29 For other data 
bases measured for other patients, the functional depen- 
dence of Z)W and D^ 2 \x) would be the same, but with 
different numerical coefficients. The order of magnitude 
of the coefficients would be the same for all the healthy 
subjects, and likewise for those with CHF (see also Wolf 
et al. 20 ). Moreover, if one analyzes different parts of the 
time series separately, one finds, (1) practically the same 
Markov time scale for different parts of the time series, 
but with some differences in the numerical values of the 
drift and diffusion coefficients, and (2) that the drift and 
diffusion coefficients for different parts of the time scries 
have the same functional forms, but with different coeffi- 
cients in equations such as (8)-(ll) . Hence, one can dis- 
tinguish the data for sleeping times from those patients 
when they are awake. 

There is yet another important difference between the 
heartbeat dynamics of the two classes of subjects: Com- 
pared with the healthy subjects, the drift and diffusion 
coefficients for the patients with CHF are very small, re- 
flecting, in some sense, their large Markov time scale tM ■ 
Large Markov times tM imply longer correlation lengths 
for the data, and it is well-known that the diffusivity in 
correlated system is smaller than those in random ones. 
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FIG. 5. Logarithmic plot of the second moment of the 
height-difference versus m, for the actual data (circles) and 
the samples regenerated by the kernel method (squares). The 
corresponding time series are plotted in Fig. 4. 



Hence, one may use the Markov time scales, and the de- 
pendence of the drift and diffusion coefficients on x, as 
well as their comparative magnitudes, for characterizing 
the dynamics of human heartbeats and their fluctuations, 
and to distinguish healthy subjects from those with CHF. 

How accurate is the reconstruction method? Shown in 
Figure 4 is a comparison between the original time se- 
ries x(n) and those reconstructed by the Langevin equa- 
tion [by, for example, using Eqs. (5), (8) and (9)] and 
the kernel method. While both methods generate se- 
ries that look similar to the original data, the kernel 
method appears to better mimic the behavior of the orig- 
inal data. To demonstrate the accuracy of Eq. (7), 
Figure 5 compares the second moment of the stochastic 
function, C2(m) = ([x(0) — x(m)] 2 ), for both the mea- 
sured and reconstructed data using the kernel method. 
The agreement between the two is excellent. However, 
it is well-known that such agreement is not sufficient for 
proving the accuracy of a reconstruction method. Hence, 
the accuracy of the reconstructed higher-order structure 
function, S„ = (|x(ti)— x(t 2 )|"), 28,29 was also checked. It 
was found that the agreement between S n for the original 
and reconstructed time series for n < 5 is excellent, while 
the difference between higher-order moments of the two 
times series, which are related to the tails of the PDF of 
the a:— increments, increases. 

The Cascade of Information from Large to Small 
Time Scales 

We would argue that if long-range, nondecaying corre- 
lations do exist in the time series, then one may not be 
able to use the above reconstruction method for analyz- 
ing them because, as is well-known, the correlations in 
a Markov process decay exponentially. Aside from the 
fact that even in such cases the method described above 
provides an unambiguous way of distinguishing healthy 
subjects from those with CHF, which we believe is more 



FIG. 6. Test of Chapman-Kolmogorov equation for 
Axi = -0.42, Axi = and Aaji = 0.42. The solid and open 
symbols represent, respectively, the directly-evaluated PDF 
and the one obtained from Eq. (1). The PDFs are shifted in 
the horizontal directions for clarity. Values of Ax are mea- 
sured in units of the standard deviation of the increments. 
The time scales T\, T2 and T3 are 10, 30, and 20, respectively. 



effective than simply analyzing the data to see what type 
of correlations may exist in the data (see below), we ar- 
gue that the nondecaying correlations do not, in fact, 
pose any limitations to the fundamental ideas and con- 
cepts of the reconstruction method described above. 

The reason is that, even if the above reconstruction 
method fails to describe long-range, nondecaying corre- 
lations in the data, one can still analyze the data based 
on the same method by invoking an important result 
recently pointed out by several groups. 20,21 ' 23-25 They 
studied the evolution of the PDF of several stochas- 
tic properties of turbulent free jets, and rough surfaces. 
They pointed out that the conditional PDF of the in- 
crements of a stochastic field, such as the increments in 
the velocity field in the turbulent flow or heights fluctu- 
ations of rough surface, satisfies the CK equation even 
if the velocity field or the height function itself contains 
long-range, nondecaying correlations. This enabled them 
to derive a FP equation for describing the systems un- 
der study Hence, one has a way of analyzing corre- 
lated stochastic time series or data in terms of the cor- 
responding FP and CK equations. We now describe the 
conditions under which such a formulation can be uti- 
lized. In such cases, one computes the Kramcrs-Moyal 
(KM) coefficients for the increments of interbeat inter- 
vals fluctuatations, Ax(r) = x(t + r) — x(t), rather than 
the time series x(t) itself. One then checks whether the 
first and second KM coefficients that represent, respec- 
tively, the drift and diffusion coefficients in a FP equa- 
tion, have well-defined and finite values, while the third- 
and fourth-order KM coefficients are small. According to 
the Pawula's theorem, 24 the KM expansion, 





0.04 



%= 5 
^=15 





■»«.» 



FIG. 7. The drift and diffusion coefficients D w (Ax) and 
(Ax), estimated from Eq. (3) for a healthy subject, follow 
linear and quadratic behavior, respectively. 



FIG. 8. The drift and diffusion coefficients D (1) (Ax) and 
x) are estimated from the Eq. (3) for typical patients 
with heart failure, and follow linear and quadratic behavior, 
respectively. 




FIG. 9. Generalized scaling analysis of a typical healthy 
subject. Structure functions S q are displayed versus S3 in 
log-log scale. 



— P(x,t\x ,t )=J2[—g^) D^(x,t)P(x,t\x ,t ) 
k=l ^ ' 



(12) 



can be truncated after the second (diffusive) term, pro- 
vided that the third- and fourth-order coefficient D^ 4 ' 
vanish, or are very small compared with the first two co- 
efficients. If so, which is often the case, then the KM 
expansion, Eq. (12), reduces to a FP evolution equation. 
In that case, a FP equation is numerically constructed by 
computing its drift and diffusion coefficients for the PDF 
P(Ax,t) which, in turn, is used to gain information on 
evolution of the shape of the PDF as a function of the 
time scale r. In essence, if the first two KM coefficients 
are found to have numerically-meaningful values (i.e., not 
very small), while the third and higher coefficients are 
small compared with the first two coefficients, the above 
reconstruction method - Eqs. (2)- (5) - are used for the 
increments of the times series, rather than the time series 
itself. 

Therefore, carrying out the same type of computations 
described above, but now for the increaments Ax(t), the 
following results are computed for the healthy subjects, 28 



larger for the healthy subjects). These arc shown in Fig- 
ures 7 and 8. Note also that these results are consistent 
with those presented earlier for the time series x(t) it- 
self, in terms of distinguishing the two classes of patients 
through their different drift and diffusion coefficients. 

The strong t— dependence of the diffusion coefficient 
£>( 2 ) for the healthy subjects indicates that the nature 
of the PDF of their increments Ax for given r, i.e., 
P(Ax, r), is intermittent, and that its shape should 
change strongly with r. However, for the patients with 
CHF the PDF is not so sensitive to the change of the 
time scale r, hence indicating that the increments' fluc- 
tuations for these patients is not intermittent. 



The Extended Self-Similarity of Interbeat Inter- 
vals in Human Subjects 



L> (1, (Aa;,r) = -(3Aa: + 0.46) x 10" 



(13) 



Lets look at another computational method for dis- 
tinguishing healthy subjects from CHF patients. The 
method is based on the concept of extended self-similarity 
L> (2) (Aa;,r) = [(l + llr" 1 ) (Ax) 2 + (5.7 + 28. 7T" 1 )] x 1(T 2 ( ESS ) f a time series. This concept is particularly useful 



whereas for the patients with CHF we obtain, 
D [1) {Ax,t) = -(1.3Aa; + 0.18) x 10~ 2 , 



(14) 



(15) 



if the time series for interbeat fluctuations (or other types 
of time series) do not, as is often the case, exhibit scaling 
over a broad interval. In such cases, the time interval in 
which the structure function of the time scries, i.e., 



S q {T) = {\ X {t + T)-x{t)\«) 



£ (2) (Ax,r) = [(5 + 5T- 1 ) (Aa 



'13 + 66T- 1 )] x 10~ 3 
(16) 



behaves as 



S q (r) 



(17) 



(18) 



Estimates of the above coefficients are less accurate for 
large values of Ax. Also computed are the average of the 
coefficients and for the entire set of the healthy 
subjects, as well as those with CHF. Moreover, is 
about jqD^ for the healthy subjects, and about ^-D' 2 ' 
for those with CHF. Therefore, the KM expansion can 
indeed be truncated beyond the second term, and the 
FP formulation is numerically justified. 

Equations (13)-(16) state that the drift coefficients for 
the healthy subjects and those with CHF have the same 
order of magnitude, whereas the diffusion coefficients for 
the given r and Aa; differ by about one order of mag- 
nitude. This points to a relatively simple way of distin- 
guishing the two classes of the subjects. Moreover, the 
T-dcpcndence of the diffusion coefficient for the healthy 
subjects is stronger than that of those with CHF (in 
the sense that the numerical coefficients of the r _1 are 



is small, in which case the existence of scale invariance in 
the data can be questioned. However, instead of rejecting 
outright the existence of scale invariance, one must first 
explore the possibility of the data being scale invariance 
via the concept of ESS. 

The ESS is a powerful tool for checking non-Gaussian 
properties of data, 31 ' 32 and has been used extensively in 
research on turbulent flows. Indeed, when analyzing the 
interbeat time series for human subjects (and other types 
of time series), one can, in addition to the r-dependence 
of the structure function, compute a generalized form of 
scaling using the ESS concept. In many cases, when the 
structure functions S q (r) are plotted against a structure 
function of a specific order, say S^ir), an extended scal- 
ing regime is found according to, 31 ' 32 



S q (r) ~ 5 3 (r)< 



(19) 



Clearly, meaningful results arc restricted to the regime 
where S3 is monotonic. For any Gaussian process the 
exponents £ g follow a simple equation, 

C 9 = \q ■ (20) 

Therefore, systematic deviation from the simple scaling 
relation, Eq. (20), should be interpreted as deviation 
from Gaussianity. An additional remarkable property of 
the ESS is that it holds rather well even in situations 
when the ordinary scaling does not exit, or cannot be 
detected due to small scaling range, which is the case for 
the data analyzed here. 

Using the ESS concept, we analyzed 33 the fluctua- 
tions in human heartbeat rates of healthy subjects and 
those with CHF, analyzed earlier with the reconstruction 
method. The results are shown in Figure 9. For a typ- 
ical healthy subject an improved scaling behavior of the 
time series is indicated by the ESS. In Figure 10 the com- 
puted scaling exponents Cq of the structure functions are 
plotted against the order q. A mono-fractal time series 
corresponds to linear dependence of Cq on q, whereas for 
a multifractal time series (, q depends nonlincarly on q. 
The constantly changing curvature of the computed C, q 
for the healthy subjects suggests multifractality of their 
corresponding time series. In contrast, Cq is essentially 
linear for the patients with CHF, indicating mono- or 
simple fractal behavior. 

It is well-known that the moments with q < 1 and 
q > 1 are related, respectively, to the frequent and rare 
events in the time series. 30,31 Thus, for the data con- 
sidered here one may also be interested in the frequent 
events in the interbeats. In Figure 11 we show the results 
for the moment q = 0.1 against a third-order structure 
function for healthy subjects and those with CHF. There 
are two interesting features in Figure 11. First, the start- 
ing point of So.i(t) versus Ss(t) is different for the data 
for healthy subjects and patients with CHF. To deter- 
mine the distance form the origin, we define, 34 

T(r = 1) = [Si ,(r = 1) + Sl(r = l)] 1 /* . (21) 

The second important feature of Figure 11 is that there is 
a well-defined r* beyond which the plot of S0.1 (t) versus 
Ss(t) is multi- valued. One can estimate r* by checking 
when S 3 (t) > S 3 (t + 1). Moreover, if we define T(r*) 
by, 

T(t*) = K 1 (t1 + S!(t*)} 1/2 , (22) 

then, there is a time scale t* such that the values of the 
third moment before and after r* are almost the same. 
Thus, the quantity r* plays a role of a local mirror on the 
time axis. In other words, locally, ^(r) for r < r* and 
t > t* has almost the same value. In Figure 11, we show 
the time scale r* and, therefore, indicate that values of 




FIG. 10. The plot of £ 9 versus q, which has a linear depen- 
dence on q for CHF subjects, but a nonlinear dependence on 
q for the healthy subjects. 

T(t = 1) and T(r*) of the interbeat fluctuations of the 
healthy subjects and patients with CHF are different. 

Table 1 presents the computed values of T x = T(t = 1) 
for both healthy subjects and those with CHF. To com- 
pute the results, we first rescaled the data sets by their 
standard deviation and, therefore, the Ti values are di- 
mensionless. The average value of T\ for the healthy 
subjects is, T(r = 1) ~ 0.5848, with its standard de- 
viations being, a ~ 0.065. The corresponding values 
for the daytime records of the patients with CHF are, 
T(r = 1) ~ 0.5077 and a ~ 0.03, respectively. Therefore, 
on average, healthy subjects possess Ti values that are 
greater than those of the patients with CHF. However, 
note that T\ for the various data sets do not have large 
enough differences to be able to distinguish unambigu- 
ously the data sets. Indeed, as Table 1 indicates, three 
of the data sets (belonging to the day- and nighttime 
records of one healthy subject, with the third one belong- 
ing to the daytime of another healthy subject) overlap. 

To develop a more definitive criterion for distinguishing 
the data for various subjects, we compute values of T(r*). 
The results are listed in Table 2. It is evident that, in 
this case, there is no overlap in the data sets. Indeed, 
the T(t*) values for the healthy subjects are larger, by a 
factor of about 3, than those for the patients with CHF, 
hence providing an unambiguous way of distinguishing 
the data sets for healthy subjects and those with CHF. 




s 3 

FIG. 11. Plot of 5*0.1 against Sa(r) for a healthy and one 
with CHF. The results indicate that the starting points are 
different for healthy subjects and those with CHF. More- 
over, for the both data set there is a well-defined r* at which 
Sz(t) >5 3 (r + l). 



TABLE I. Values of Ti for the healthy subjects (day time 
and night time ) and those with CHF. 



Healthy 


CHF 


0.658 


0.557 


0.672 


0.565 


0.614 


0.539 


0.605 


0.526 


0.583 


0.512 


0.581 


0.493 


0.576 


0.492 


0.558 


0.481 


0.494 


0.469 


0.480 


0.443 







Comparison with other Methods 

Stanley and colleagues 5 ' 13,19 ' 26 ' 27 ' 38 and other 35-37 an- 
alyzed the type of data that we consider in this review 
by different methods. Their analyses indicate that there 
may be long-range correlations in the data, characterized 
by self-affine fractal distributions, such as the fractional 
Brownian motion, the power spectrum of which is given 
by 

S(f) ~ r (2H+1) , (23) 

where H is the Hurst exponent that characterizes the 
type of the correlations in the data. Thus, healthy sub- 
jects are distinguished from those with CHF in terms of 
the type of correlations that might exist in the data: neg- 
ative or antipersistent correlations in the increments for 
H < 1/2, as opposed to positive or persistent correlations 



TABLE II. Values of the T(r*) for the healthy subjects 
(day time and night time) and for those with CHF. 



Healthy 


CHF 


3.08 


0.741 


2.68 


0.714 


2.34 


0.685 


1.92 


0.681 


1.86 


0.675 


1.42 


0.632 


1.40 


0.573 


1.22 


0.728 


1.20 


0.552 


1.15 


0.465 







for H > 1/2, and Brownian motion for H = 1/2. While 
this is an important and interesting result, it can also be 
ambiguous and not very precise. For example, suppose 
that the analysis of two time series yields two values of 
H, one slightly larger and the second one slightly smaller 
than 1/2. Then, 2. Its thus difficult to state with confi- 
dence that the two times series are really distinct. 

The reconstruction method described earlier analyzes 
the data in terms of the Markov processes properties. As 
a result, it distinguishes the data for healthy subjects 
from those with CHF in terms of the differences between 
an FP equations drift and diffusion coefficients. Such dif- 
ferences are typically very significant and, therefore, pro- 
vide, in our view, an unambiguous way of understanding 
the differences between the two groups of the subjects, 
including for those series for which the Hurst exponents 
arc only slightly different. In addition, the computational 
approach described in this review provides an unambigu- 
ous way of reconstructing the data, hence providing a 
means of predicting the behavior of the data over periods 
of time that are on the order of their Markov time scales. 

Although it remains to be tested, we believe that, to- 
gether, all the computational methods that have been 
described in this review are more sensitive to small differ- 
ences between the data for the two groups of the subjects 
and, therefore, might eventually provide a diagnostic tool 
for early detection of CHF in patients. 

Finally, the computational approaches described in 
this review arc quite general, and may be used for analyz- 
ing times series that represent the dynamics of completely 
unrelated phenomena. For example, we have used the 
concepts of Markov processes and extended self-similarity 
to develop 34 a method for providing short-term alerts for 
moderate and large earthquakes, as well as making pre- 
dictions for the price of oil. 39 
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